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Abstract. Recently, Narayan & Raymond (1999) pro- 
posed that the thermal emission lines from the hot plasma 
in advection dominated accretion flows (ADAFs) are po- 
tentially observable with the next generation of X-ray ob- 
servatories, with which the physical properties of some 
X-ray sources can be probed. In ADAFs, the tempera- 
ture of the ion is so high that the thermal broadening of 
the line is important. We calculate the profiles of ther- 
mal line emission from ADAFs, in which both the ther- 
mal and Doppler broadening have been considered. It is 
found that the double-peaked profiles are present for high 
inclination angles between the axis of disk and the line 
of sight. The double-peaked profiles are smeared in low 
inclination cases, and completely disappear while the in- 
clination angle is less than 15°, where the thermal and 
turbulent broadening dominated on the line profiles. We 
also note that the thermal line profile is affected by the lo- 
cation of the transition radius of ADAF. The self-similar 
height-integrated disk structure and the emissivity with 
power-law dependence of radius are adopted in our calcu- 
lations. The results obtained in this work can be used as a 
diagnosis on the future X-ray observations of the thermal 
lines. Some important physical quantities of ADAFs could 
be inferred from future thermal line observations. 
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1. Introduction 

Emission line profiles from accretion disks have been inves- 
tigated by many researchers from both observational and 
theoretical sides. In the past decade, the X-ray line emis- 
sion has been extensively studied, mainly stimulated by 
the progresses on X-ray observations. The profiles of emis- 
sion line from accretion disks around a Schwarzschild black 
hole (Fabian et al., 1989; Matt et al., 1993), or around a 
Kerr black hole (Laor 1991) are well calculated. The line 
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spectra from ADAFs have been calculated by Fukue & 
Ohna (1997) without considering the thermal broadening 
of line profile in the rest frame of the emitter. The fluores- 
cent line emission from cool gas irradiated by hot corona 
has been studied in some detail (Mushotzky et al., 1993; 
Tanaka et al.,1995), which is the focus in most X-ray line 
formation mechanism investigations. The fluorescent line 
can be formed in the disk region very close to the black 
hole. However, the thermal X-ray line emission from hot, 
optically thin gas around black holes has not been dis- 
cussed very much. More recently, Narayan & Raymond 
(1999) present the calculations on thermal line formation 
in hot ADAFs, and propose that such X-ray emission lines 
are potentially observable with the next generation of X- 
ray observatories. The physical properties of some X-ray 
sources and the constraints on the accretion flows in these 
sources may be provided by future observatories. 

In a two-temperature ADAF, while the ions receive 
most of the viscous energy and are nearly virial, the elec- 
trons are heated by several processes, such as, Coulomb 
coupling with the ions, compression, and direct viscous 
heating, and cooled by a variety of radiation processes 
(Narayan & Yi, 1995b). The temperature of the electron 
is in general lower than that of the ions due to insufficient 
interaction between electrons and ions(see Narayan et al. 
1998 for a review on ADAFs). The thermal X-ray line 
emission is supposed to be from the disk region where the 
electron temperature less than 10^ K, which corresponds 
to the disk region r > lOOrg (Narayan & Raymond 1999). 
The profiles of the fluorescent line are mainly determined 
by the Doppler broadening caused by the bulk motion of 
the flows, and the thermal broadening can be neglected in 
the profile calculations. Unlike the case of the fiuorescent 
line, the thermal X-ray lines are supposed to be from the 
disk region r > lOOrg with the ion temperature around 
10^" K. In this case, the thermal broadening of the line 
profiles is important compared with the Doppler broaden- 
ing, since the bulk velocity of the ADAFs is only a fraction 
of the Kepler velocity, and the Kepler velocity at lOOrg is 
only one-fifth of that at the inner edge of the disk around 
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a Schwarzschild black hole. In this work, we present calcu- 
lations on the thermal line profiles from ADAFs. In Sects. 
2-4, we give the formalism of the problem. The last section 
contains the results and a discussion. 

2. Thermal broadening of the line profiles 

We first calculate the thermal broadening of line emis- 
sion from relativistic ions. Using the common relativistic 
transformations and assuming the motion of thermal ions 
to be isotropic, we write the observed flux density of the 
line from the ions moving at /3i = vij c as 
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where eiine is the line emission power in the rest frame of 
the ion, v^ Vp are the frequency measured by the observer 
and the line frequency measured in the emitter's frame 
respectively, 7i is the Lorentz factor of the ion's thermal 
motion, and /i — cos Q (Q is the angle between the motion 
of the ion and the line of sight). Using the relation Vf^ — 
v^i(\ — Pifi), we can rewrite Eq. (1) as 
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where 0ion is the thermal temperature in units of ion rest 
mass, and K2 is the second-order modified Bessel function. 
The ion distribution function is normalized as 



N{^i)d-i, = 1. 



We can therefore obtain the broadening profile of the ther- 
mal line emission as 
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For the given ion temperature Qjon, we can integrate 
Eq. (4) and obtain the line profile: 
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Besides the thermal broadening, the turbulent velocity 
of eddies in ADAFs may also affect the line profile. We 
assume the turbulence to be isotropic and adopt a mean 
turbulent velocity vt = (3tc. The thermal line profile can 
be calculated by 

Lilt, Bion, t^) 



4A7t^2(l/eion)^^| 



Bio 



2x2 



exp 



1 



■^^ion-^ 



dx. 



(6) 



where 
and 
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3. Structure of self-similar ADAF 

We consider the self-similar disk structure given by 
Narayan & Yi (1994). In this solution, the azimuthal ro- 
tational velocity w^ and radial velocity Vr are given as a 
function of r as follows: 
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It is assumed that magnetic fields contribute a constant 
fraction (1 — /?) of the total pressure in ADAFs. The ap- 
propriate relation between the ratio of the specific heats 
7 and P is given by Esin (1997): 7 = (8 - 3/3)/(6 - 3/3). 
The disk structure is then described by three parameters: 
the ratio of the magnetic pressure /3, the viscosity a, and 
the fraction / of viscously dissipated energy which is ad- 
vected. 

For simplicity, we further assume that the emissivity j 
of the unit area of disk surface has power-law dependence 
of radius r, 



j{r) ex r 



(9) 



The radial velocity Vr and azimuthal velocity v^ of 
the flow given by Eqs. (7) and (8) will be used in the 
calculations of line profile. 
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4. Profiles of thermal lines from ADAFs 

Assuming the line of sight to be inclined at an angle i to 
the disk axis, the angle S between the line of sight and the 
motion of gas in the disk is 



cos 6 = sin i cos 



(10) 



where tan0„ = Vr/v,p. Using the relativistic transforma- 
tions, the line profile for a ring in the disk between r and 
r + dr is available. 
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(12) Fig. 1. The velocity of the accretion flow as functions 
of /3 for different values of viscosity: a =0.1 (dashed), 



The electron temperature Te has been neglected in Eq. 0.2(solid), 0.3(dotted), 0.5(dash-doted), respectively. 
(12) (see Eq. (2.16) in Narayan & Yi 1995b), since T^ 
is about an order of magnitude lower than Ti in two- 
temperature ADAFs. In this paper, we focus our calcu- 



lations mainly on iron, i.e., Qi, 
between v' and z^„h, is 
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where the factor (1 — 1/r) ^/^ represents the gravitational 
redshift. From Eqs. (10) and (13), we can obtain 
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Finally, we can calculate the line profile as 

•Tout 

Fy{vohs) = / rjdr COS i 



Fig. 2. The thermal iron line profiles with parameters: 
a = 0.1, 7 = 1.444 (P = 0.5), and i = 0°, 60°, respectively. 
The dotted lines indicate the profiles without considering 
turbulent motions in ADAFs. 
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The inner and outer radii of the line emission region are 

assumed to be ri„ — lOOr^ and rout — lOOOr^. The typical 

T/2^^ ) turbulent velocity of eddies in ADAFs is approximated as 

/ ,^ \i/2 
vt ~ acs, Cs = i]]-^] , where ^i = 1.23 is adopted, 

same as Narayan & Yi (1995b). Now, using Eqs. (6) and 
(15), we can calculate the profiles of the thermal line from 
ADAFs. 



5. Results and discussion 

We depict the total bulk velocity of the flow as functions 
of the disk parameters /3 and a in Fig. 1. It is found that 
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Fig. 3. Same as Fig. 2, but a = 0.2. 




Fig. 4. Same as Fig. 2, but a = 0.5. 



the velocity is mainly determined by the value of /3. The 
viscosity a affects the velocity very little, especially in the 
range of typical values of a around 0.2 ~ 0.3 (Narayan et 
al. 1998). One may therefore draw a conclusion that the 
line profile is mainly determined by the value of magnetic 
pressure ratio j3 that determines the ratio of specific heats 
7 here. We note that Quataert & Narayan (1999) give a 
different relation between the ratio of specific heats 7 and 
plasma /3, in which the particle heating is explicitly in- 
cluded. The plasma /3 is a function of the radius r, and 
the structure of ADAFs is therefore sensitive to the outer 
boundary condition. In this work we limit our calculations 
to Narayan & Yi's (1994) self-similar ADAF model in or- 
der to get the basic features of thermal line profiles. The 
present thermal line profile calculations can be easily ex- 
tended for different ADAF models. The observed thermal 
line profiles would be used as tests on different ADAF 



Fig. 5. The profiles of the thermal iron line emission from 
accretion flows for model A: a = 0.2, / = 1, 7 = 1.444 
(/3 — 0.5), with different inclination angles: i = 0° (solid), 
15° (dashed), 30° (dotted), 60° (dash-dotted). 




0.96 



Fig. 6. The profiles of the thermal iron line emission from 
accretion flows for model B: a = 0.2, / = 1, 7 = 1.602 
(/? = 0.9), with different inclination angles: i = O°(solid), 
30° (dashed), 45° (dotted), 60° (dash-dotted). 



models, and some physical quantities of ADAFs can be 
inferred from X-ray observations on thermal line profile 
and continuum emission. 

In Figs. 2-4, we compare the importance of turbulent 
motions in the formation of the line profiles for different 
values of viscosity a. For small a, the turbulent velocity 
of eddies is small, and the line profiles arc hardly affected 
by turbulence, while the profiles are obviously broadened 
in high a cases (see Fig. 4). 

We adopt two sets of the disk parameters: A. a = 0.2, 
/ = 1, /3 == 0.5; B. a = 0.2, / = 1, /3 = 0.9. The ratio of 
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Fig. 7. The profiles of thermal iron line emission from Fig. 9. Same as Fig. 5, but for thermal oxygen lines, 
accretion flows for model A with i — 60° for different 
emissivity law: h ~ (solid), 1 (dashed), 2 (dotted). 
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Fig. 8. Same as Fig. 7, but for model B. 



Fig. 10. The profiles of thermal oxygen line for model 
B, with different inclination angles: i — O°(solid), 
45° (dashed), 60° (dotted), 75° (dash-dotted). 



specific heats is 7 ~ 1.444 for j3 — 0.5, while 7 ~ 1.602 for 
(3 — 0.9. The line profiles for models A and B are plotted in 
Figs. 5 and 6, respectively. The exponent of the emissivity 
law is assumed to be 6 = 1 in these calculations. We find 
that the double-peaked profiles appear for high inclina- 
tion angles and are smeared by the thermal and turbulent 
broadening in low inclination cases. The double-peaked 
structure finally disappears when the line of sight is close 
to the axis of the disk. The bulk velocity of accretion flows 
is large in the low-/? case, and one can then find in Figs. 
5 and 6 that the wings of line profiles for model A extend 
over a larger range of frequency than that for model B. 
The double-peaked structure disappears completely when 
the inclination angle is about 15° for model A, while an 



obvious plateau is present at the inclined angle around 45° 
for model B (see Fig. 6). 

In Figs. 7 and 8, the line profiles corresponding to the 
different emissivity law 6 = 0, 1, 2 are plotted for mod- 
els A and B respectively. The line emission from the outer 
region of the disk dominates in the case of 6 = and there- 
fore has a narrower line profile, while a broader line profile 
is present in the case of 6 = 2, since the line emission is 
mainly from the inner region of the disk in this case where 
both the bulk and thermal velocities of the flow are high. 

The line profile calculations in this work are mainly 
performed for iron, which would be most probably ob- 
served by X-ray observatories. We also plot the thermal 
oxygen line profiles in Figs. 9 and 10 for models A and B, 
respectively. It is found that the double-peaked structures 
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Fig. 11. The line profile dependence of the outer ra- 
dius for model A: rout = SOOr^ (dashed). lOOOrg(sohd), 
5000rg (dotted). The inchnation angle is i = 45°. 




Fig. 12. Same as Fig. 11, but for model B. 

are smeared even for high inchnation cases. For model B 
(7 — 1.602), the bulk velocity of ADAF is relatively low, 
the double-peaked structure does not appear for any in- 
clination angle, and only a plateau is present for the high 
inclination angle (see Fig. 10). The thermal lines for light 
ions have systematically broad profiles due to the fact that 
the lighter ion has higher thermal motion velocity for the 
same ion temperature. The expected difference between 
line profiles for different ions can be tested by future X- 
ray observations. 

It is proposed that a transition radius rtr exists in the 
ADAF model. In the region of r < rtr, the accretion flow 
is a two-temperature ADAF, while in the region beyond 
the radius rtr, the accretion occurs partially as a thin disk 
(Narayan, McClintock & Yi 1996). The profiles of thermal 



line will be obviously affected by the transition radius rtr, 
which may be related to the accretion rate m (Narayan 
& Yi 1995b). In present line profile calculations, we have 
fixed the outer radius of line emission region as lOOOrg. 
We depict the line profile dependence of the outer radius 
in Figs. 11 and 12. It is found that the line profile becomes 
narrow and a small difference between two peaks is present 
for a large outer radius. The outer region of the ADAF has 
relatively lower bulk velocity, and the emission from this 
region contributes mainly on the central part of the profile, 
which leads to a narrow line profile for a large outer radius. 
In the case of the ADAF with a small outer radius, the 
Doppler beaming effect caused by bulk motion of the flow 
in the whole line emission region is strong, and therefore 
the line profile is broad. 

The line profiles are calculated on the assumption of 
the thermal distribution for ions (but also see Mahadevan 
& Quataert 1997). The thermal and turbulent broadening 
of the line is important, but the Doppler broadening by 
the motion of accretion flow cannot be neglected either. 
The profiles for thermal lines from ADAFs show that the 
blue peak is higher than the red one and line profiles ex- 
tends more in the red wing, which is similar to previous 
calculations on the line profiles from the standard accre- 
tion disk. The difference between two peaks of the line 
from ADAFs is smaller than that of standard accretion 
disk cases, and the peaks are finally smeared by the ther- 
mal and turbulent broadening when the inclination is low. 
We note that the profiles of the neighbouring X-ray lines 
calculated by Narayan & Raymond (1999), for example, 
FexxvAl.855 and FexxviAl.780, may partly overlap in 
high inclination cases. The present calculations on ther- 
mal line profiles can be used as a diagnosis on the future 
X-ray observations of the thermal lines. 

We have not considered the Comptonization in ADAFs 
in present calculations, which may affect the thermal 
line profiles to some extent. It can be included in future 
line profile calculations, perhaps by numerical approaches 
(Dumont et al. 2000). 

Our line profile calculations are based on the self- 
similar height-integrated ADAF models. The vertical 
structures of the velocity, density and ion temperature 
have not been taken into account. The velocity field in 
ADAFs may depend on the vertical height (Narayan & Yi 
1995a; Igumenshchev & Abramowicz 1999). Narayan & 
Yi (1995a) studied the structure of ADAFs in rO plane by 
numerically solving the self-similar equations. They com- 
pared their exact self-similar solutions with those obtained 
using the system of height-integrated equations and found 
that various dynamical variables such as the radial ve- 
locity, angular velocity and sound speed estimated from 
height-integrated approximate solutions agree very well 
with the corresponding spherically averaged quantities in 
their exact solutions. Their exact solutions also show the 
angular velocity almost not varying with the angle at 
given radius r for small e'. The present calculations can 
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therefore be good approximations of the reality, provided 
the emissivity is interpreted as that from the shell at given 
radius r. We believe that the main features of line pro- 
files have not been affected much by the hmitation of 
one-dimensional velocity structure adopted in the calcula- 
tions. The further calculations considering rO distribution 
of variables Vr, v^, p and ion temperature Ti are necessary 
for detailed modelling on line profiles. 
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